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Abstract 

We show that the presence of finite-size monopoles can lead to 
a number of interesting physical processes involving quantum entan- 
glement of charges. Taking as a model the classical solution of the 
N = 2 SU (2) Yang- Mills theory, we study interaction between dyons 
and scalar particles in the adjoint and fundamental representation. 
We find that there are bound states of scalars and dyons, which, re- 
markably, are always an entangled configuration of the form = 
|dyon + ) |scalar_) ± |dyon_)|scalar + ). We determine the energy levels 
and the wave functions and also discuss their stability. 
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1 Introduction 



The presence of a 't Hooft-Polyakov monopole [|l], § gives rise to remark- 
able effects such as baryon number violation H, Hfl. It is also believed that 
there must exist processes where a charged particle transfers charge to the 
monopole, converting it into a dyon (see e.g. |5], ||, [7|, U), or more generally, 
that there are charge transfer processes in the interaction of particles with 
dyons. These effects are believed to be quite generic, peculiar to the nature 
of the monopole field, but otherwise model independent. In particular, they 
are present in supersymmetric models. 

Here we will consider specifically the Prasad-Sommerfield dyon solution of 
N = 2 supersymmetric SU(2) Yang-Mills theory (for reviews see |I| |TT|1 ) 
and its interaction with scalar particles of a SU(2) multiplet in the adjoint 
and in the fundamental representation. 

In general, the equations of the scalar field in the background of the dyon 
couple the different members of the SU(2) multiplet. By defining the electric 
charge in terms of the asymptotic states where the Higgs field is constant, 
far from the dyon core, this fact is seen as a nonconservation of the charge 
of the scalar field. Since it is possible to define the total electric charge by 
the flux at infinity of the total electric field (that is, including the electric 
field of the dyon solution) and show that it is conserved, this implies that 
in the process there is a transfer of charge from the particle to the dyon 
core. While it is difficult to incorporate exactly the back reaction of the 
scalar field on the dyon solution, it is possible to take into account the main 
effect of the total charge conservation by representing the charge degrees of 
freedom of the dyon by means of a quantum rotator formalism (1)12, 11] , see 
in particular ref.|J). The outcome is that in general the quantum states of 
the particle-dyon system are of the form of quantum entanglement, i.e. a 
linear superposition of particle-dyon states in which the total charge is fixed, 
but the particle and the dyon appear with various charge assignements. 

It is interesting to compare with the string theory description of the 
N = 2 monopole in terms of D branes. The N = 2 SU(2) monopole can be 
geometrically described as a tube (representing a D string) connecting two 



parallel D3 branes |13[ . The vacuum expectation value of the Higgs field is 
proportional to the distance between the D3 branes. This separation breaks 
the SU(2) gauge group to U(l). Consider a small open string representing a 
neutral particle on one of the two D3-branes: one of the string end points, say 
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the negatively charged one, can fall into the monopole and pass through it to 
the other brane, whereas the other endpoint may remain on the first brane. 
An observer on the first brane would see a positively charged particle and a 
negative charge dyon, due to the flux lines of the electric field dragged by the 
particle which failed into the monopole. The string would go from the positive 
particle to the dyon (providing a physical realization of the mathematical 
string introduced in ref . || for preserving gauge invariance) and pass through 
the tube up to the negative end point on the other brane. Clearly, it is 
equally probable that the positively charged endpoint of the original open 
string falls into the monopole, which would lead to a system of a positive 
charge dyon and a negative charge particle on the first brane. Of course 
this geometrical picture is unable to take into account the quantum effect 
of entanglement, whereby the resulting state can be a superposition of the 
two possible realities: a dyon (+) and a particle (— ), or a dyon (— ) and a 
particle (+). 

In this paper we consider the possible bound states of a dyon and a 
charged scalar particle. We first perform a general harmonic analysis in the 
case of the finite-size monopole and, after diagonalizing the angular momen- 
tum, we get coupled radial equations, where the off-diagonal terms - repre- 
senting the coupling of different charges - vanish exponentially outside the 
dyon core. As a result, we find that the possible bound states are always 
of the entangled form. In other words, a bound state just of the form, say, 
|dyon_) |scalar + ) (a plus charged scalar particle and a minus charged dyon) 
can never occur. Rather, the possible bound states are a linear combina- 
tion of the above state with the state of opposite charge, |dyon + ) |scalar_). 
At large distances, the system is similar to a hydrogen atom, with a parti- 
cle orbiting far away from the core. We will find that the energy spectrum 
and wave functions can indeed be approximated by the same formulas of 
the hydrogen atom, with a proper definition of the effective charge and ef- 
fective angular momentum. The most relevant effect of the core is thus the 
production of quantum entanglement. 

Bound states of particles with dyons have been studied in the past (see 
e.g. [|14], 0), mainly for fermions and in the point-like core limit, except 



for the special zero energy state found by Jackiw and Rebbi ||16|| . However, 
the peculiar fact that bound states to dyons are always quantum entangled 
with the system of opposite charges does not emerge in this limit. In ref. 
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I7|, 113] a more refined analysis was performed, by keeping into account the 



core interior in some stepwise approximation, but focusing on the case of the 
monopole or for a fixed dyon, and mainly for investigations on the Jackiw- 
Rebbi phenomenon; dynamical effects of charge transfer and related collective 
dyon modes were not taken into account. 

In order to observe the quantum entanglement effect, the key points are: 
a) the interactions in the monopole interior (in particular, the effect would 
not appear in a large distance approximation neglecting what happens inside 
the core); b) the account of the dyon degrees of freedom and of the collec- 
tive mode which is responsible for charge conservation. Here we report on 
this phenomenon, not only by controlled analytic approximations but also 
providing the exact numerical solutions to the bound state problem. 

Our computation is done for the case of globally neutral systems of scalars 
and dyons, for scalars in the adjoint and fundamental representation of 
SU(2), for the case of scalars belonging to the N = 2 supermultiplet and 
also more in general. We also discuss briefly the stability of the bound states 
against perturbations not included in our computation, in particular, the 
radiation of e.m. quanta. A summary of the results is reported in the last 
Section 5. 



2 Dyons in N = 2 577(2) super Yang-Mills 
2.1 Monopole and dyon solutions 

Let a a , a = 1, 2, 3 be the Pauli matrices, satisfying [a a , a b ] = 2ie abc a c . The 
background describing the dyon of N = 2 SU (2) super Yang-Mills theory is 
given by @, ^ 

A 1 = A\o h = -o h e hl3 — 2 fix) , <$> b a b = —aG(x) , (1) 
er z r 

T (J 

A = -^—aG(x) sin 6 , x = ear cos B , (2) 

f{x) = - K(x)) , K(x) = — ^ , G(x) = cotha; - - , 

smn(a;j x 
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where a represents the expectation value of a scalar field, and (ea cos 0) 1 
determines the dyon size. The behavior at large x and small x is 



K(x) = 0(e~ x ) , G(x) = 1 - - + 0(e~ x ) , x > 1 , 

Ob 



K(x) 



x 



1 - — + 0(x 4 ) , = | + 0(a; 3 ) , x < 1 

o 6 



The dyon charge and mass are given by 



An 



q = g tan = tan , q = ne 

e 



M = aJg 2 + n 2 e 2 



ag 



cos 



For e 2 C 1, this becomes 



M = \ag\(l + \a 2 ) , a < 1 , 



a = tan = — 



eg 

47T 



(3) 
(4) 



(5) 

(6) 
(7) 



In particular, the mass difference between a dyon and a monopole is given 
by 



AE = M dyon - M, 



mon — 2 



ago; | 



(8) 



The Coulomb potential has a constant piece at infinity, related to the energy 
associated with the dyon electric field. For the analysis in Section 3, it is 
convenient to write it as 



A o = !jZ Vq y = -^G(x)cos0 . 



(9) 



where q represents the dyon charge operator. 



2.2 Charge conservation 

Let $ = $5(j 6 be the Higgs field which has a constant expectation value at 
infinity, and let ^ n stand for the other scalar fields. We choose the gauge 
A = 0. 
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The Lagrangian is invariant under the transformation 

5$ = 0, 5^ n = z[e$, , <^z = e$$ + i[e$,A,] • 

Consider now e = e(t) as an arbitrary function of time. The Lagrangian is 
no longer invariant but the variation of the Action with respect to e must be 
zero because of the equation of motions. We thus find a conserved quantity 
Q: 

Q= f d 3 r Ti(D,$^+i^$[$ B) ^ B ]) . (10) 

One can indeed independently check that d t Q = by using the field equa- 
tions. Note that Q is given by an integral of the sum of two densities, the 
"charge density of the scalar fields" \& n , that is iTr J2 n &l^n, d t ^ n ] plus the 
"charge density of the dyon", that is Tr(Di&d t Ai). A similar proof for charge 
conservation is given in || for the case of SU(2) Yang-Mills theory with an 
isodoublet of Dirac fermions. 
Using the field equation 

5M = iE[*n,^n]+i[$ l ^] , 

n 

the charge can be written as 

Q = J dETr[<f>d t A] , 

where S is a closed surface at infinity. 

Thus Q is interpreted as the total charge (apart from a constant factor, 
equal to the expectation value of the Higgs field at infinity) since it is mea- 
sured by the total electric flux in the spontaneously broken theory. Quantum 
mechanically, if the initial state is an eigenstate of Q, also the final state must 
be an eigenstate of Q, since Q commutes with the Hamiltonian. 

3 Scalar particles in the adjoint representa- 
tion 
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3.1 Covariant Equations 

Let us consider a scalar particle in the dyon background. The scalar particle 
is a quantum of the matrix valued scalar field we take an SU(2) triplet 
^> = ^aCr". The equation of motion is given by 

^■* = ie 2 [$,[$,*]] , (11) 

where D\ ■ * = D„ o (D„ o ), D„ o = - ie[A„, % 

We look for solutions of eq. ( |i~l~l) of the form: *&(r, t) = e~ lEt ijj(f), and 

ip{r) is interpreted as the stationary wave function of the scalar particle. 
It is convenient to choose the following basis for SU(2) Lie algebra: 

a = UiOi , (3 = ViOi , 7 = n i°i > ( 12 ) 

Ui = 5 3i - n 3 ni , Vi = e 3ji nj , m = — . (13) 

Note that — , — , n, with n% = 1 — n\ represent an ortho normal frame for 

TLX TIT 1 

vectors in R 3 , u.v = u.n = n.v = 0, n 2 = 1, u 2 = v 2 = n\. It follows that 
Tr a$ = Tr d 7 = Tr 7( 5 = 0, ±Tr a 2 = |Tr (5 2 = n 2 T} ±Tr 7 2 = 1. 

The matrices d,/3, 7 obey the commutation relations: 

[a, 7] = 20 , [(3,y) = -2i&, [a, ft = -2m| 7 , 

and 

[d±, j] = ±2a± , a± = a ± i/3 . 
The matrix valued wave function if)(r) can be decomposed as 

ip{r) = F+(r) a- + F_(f) a+ + F (f) 7 . (14) 

The component F (r) multiplying 7 represents a neutral component, whereas 
the components F± of a =F i/3 represent charge e and charge — e components, 
respectively. This is clear by a gauge transformation so that 

U~ 1 jU = <t 3 , U- 1 (a±t(3)U = -2sinfl e ±ilp a T , (15) 

where 

U(n) = cos- +?sm- — 7— -p . (16) 
6 



In this gauge, the charge operator is Q = | 03, acting on a state \I/ by 

q. * = [&*]. 

Now we evaluate .D 2 ■ ■0. With the definitions (1.4), and further defining 
for short Q = d 2 - 4(1 + K 2 ) - 2 -^-u j d h we find: 

D 2 -(iF (r) )=^(d 2 - 2 -^-)F + a (^^.)f + /3 (^^)f . (17) 



D 2 ■\{a±if3)F T ] = (a±i(3)(tt±-^iv j d j )F T 



r n 



T 



+ l 2 ^{T^v ] ^ ] -u 3 ^ J + 2 ^y T . (18) 
3.2 Harmonic analysis 

To solve for the angular dependence, it is convenient to use spherical coor- 
dinates r, 9, ip. We have: 

Ujdj = — sin 9 dg , ivjdj = i-d v . (19) 
Asymptotically at r — > 00, we have K — > and in this case there is decoupling 



of charges (+,—,0) in fll7| ) and ([Tq), as expected. Consider the monopole 



case = 0. Charged particles have mass m, with m = ea. Define 



F±{f) = —,h ± (r) . (20) 
9 



We find that for r 3> m 1 eqs. ([□]), ( [18] ) give 



—d r r 2 d r + 1 <9g (sinfldg) + - 1 2 (<9 y =f i(l - cos6»)) 2 
r 2 r 2 sm0 ran ^ 



/ - (m- -) 2 )// . (/•') = . 



This is the covariant equation for a scalar particle of charge =Fe moving in a 
£7(1) point-like Dirac monopole background. In addition, there is a Coulomb- 
like potential due to the long range tail of the Higgs field. 
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Let us now consider the complete equations in all the space. In the 
Appendix we report the harmonic analysis in full detail. The angular depen- 
dence is solved by setting 



Z = y/l(l + 1) 



where Yi m (6,(p) are the standard spherical harmonics, and can be ex- 
pressed in terms of them: 



e Tiip 

ZL = —^(smOde =F id v )Y lm . (21) 
sin u 

Note that the angular expansion of the charged components h± begins with 
1 = 1. 

Using the results of the Appendix, we get the following system of coupled 
differential equations for 4>o(r), 4>±{r) : 

- ^r'dr + + ^Kl + l))0o + ^(0+ + M = E 2 <Po , (22) 



^r 2 ^ + ^(_i + K 2) + J_/(/ + i)] 0+ + 



= [(£ - eVq) 2 - e 2 a 2 G 2 ]<f) + , (23) 

[ - ^d r r 2 d r + £(-1 + if 2 ) + i/(/ + 1)]0_ + 

= [(E + eVq) 2 - e 2 a 2 G 2 ](f)- , (24) 

where V was defined in eq. (|9]). 

3.3 Collective coordinates and equations of motion 

When a process of charge transfer occurs, there is an important back re- 
action in the dyon background. This effect can be taken into account by 
incorporating into the dynamics the collective coordinate x associated to 
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global £7(1) gauge transformations. This "quantum rotator" degree of free- 
dom ensures charge conservation. The dyon charge operator is written as 
q = gtan© = —ied x . 

Consider a system consisting of a scalar particle and a dyon with total 
charge equal to zero. The wave function can be written in the form 



/ 6+ e~ ix ' 



(25) 



V 



In general, this represents a mixture of monopole and neutral scalar, with 
dyons and scalars of charges (H — ) and ( — h). When the system makes a 
transition from a monopole to a dyon, the field component A Q is turned on. 
In addition, the energy of the scalar is reduced into an amount equivalent 
to the mass difference between the dyon and the monopole. In what follows 

2 

we will assume that a = f- is small. If a is not small, other back reaction 

4-7T ' 

effects become important and this semiclassical analysis is not applicable (for 
a < 1, one may approximate sinB = tan0, cosG = 1 + 0(e 4 ) = 1). 
The equations of motion (p2|) -(p4|) can be written as 



where 



Hip = 0, (26) 



H = ( P + (E + msGid x f - m 2 G 2 ) l 3x3 + %M 



M 



I Ke~*x 

Ke ix Ke~ ix I . .s = sin (-) 
V Ke ix 



1 1 IK 2 1 

P = — 2 d r r 2 d r + -/(/ + 1) + — + -.(1 + K 2 )d\ . 

Then the equations of motion for the neutral system take the form 



J_ 5r r 2 a r + ^- + ii(J + l))0o + ^(0+ + <P-) = E 2 <p , (27) 



[ - ^d r r 2 d r + £(-1 + K 2 ) + ±1(1 + 1)]0 + + ^a0 o 

= [{E - msGf - m 2 G 2 ]<f) + , (28) 
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[ - ^d r r 2 d r + £(-1 + K 2 ) + ^1(1 + 1)]0_ + 

= [(£ - msGf - m 2 G 2 ](j)- . (29) 

The incorporation of the collective coordinate \ produces the flip of sign of 
the dyon charge in the equations for + and 0_, required by charge con- 
servation. As a result, the Coulomb potential of the dyon-scalar is always 
attractive for both + and 0_. 

In terms of the variables a = 0+ — 0- and S = 0+ + 0_, we have a 
decoupled equation 

[ - ^ r r 2 9 r + £(-1 + K 2 ) + + l)]0 a 

= [{E - msGf - m 2 G 2 ]0 a , (30) 

and a coupled system of ordinary differential equations for S and 0o, 

[ - ^d r r 2 d r + £(-1 + K 2 ) + 4,i(Z + 1)]0 S + ^0o 

= [(£ - msGf - m 2 G 2 ]0 s , (31) 

( - ~ 2 d r r 2 d r + % + + 1))0 O + ^0 S = £ 2 o . (32) 

In the Section 3.4, we will see that eq. (|3H| ) describes bound states. The 
wave function for such bound state is of the form = |dyon + ) | scalar _) — 
|dyon_) |scalar + ), i.e. there is quantum entanglement of charges with a dyon 
state. Eqs. fl31]) , (|32|) can be used to describe the dynamics of a scattering 
process in which a neutral scalar particle scatters off a monopole and gives 
rise to an outgoing state which is a mixture of a neutral scalar and a +/— 
charged scalar entangled with a — /+ charged dyon. 



3.4 Bound states 

Here we study the solutions of the equation (|30"D for a . In terms of the radial 
coordinate x = (mcosG) r, it takes the form 

[ - ±zd x x 2 d x + £(-1 + K 2 (x)) + ±1(1 + l)]0 a 

= [£ 2 - 2£aG(x) - G 2 (x)]0 a , (33) 

E ^ e 2 

£ = , a = tan 6 = — . 

m cos 47Y 
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This is a Schrodinger-type equation with a potential which is an attractive 
Coulomb potential at large distances, where G = 1 — - (see eq. (0) ), and 
therefore it may admit bound state solutions. 

A closed analytic solution of this differential equation is not known. We 
will solve it by using two independent methods: a) Analytic method; b) 
Numerical method. 

Let us begin with the analytic method. We approximate the potential 
by a simpler function, replacing K 2 and G by Kq, Go defined as follows (cf. 
eqs. (D,(|) ): 

K 2 (x) = , G (x) = 1 - - , x>l , 

X 

.2 

Kl{x) = 1 - y , G {x) = , X < 1 . 
Exterior solutions: For x > 1 the equation is 

[ - - 2 d x x 2 d x + - - (a£ + 1) + k 2 ]<p a = (34) 

with 



fe = Vl + 2a£ - £ 2 . (35) 

Note that the Coulomb potential in eq. (0) has two contributions, one 
coming from the interaction with the Higgs field, the other from A Q . For 
a C 1, the Coulomb potential due to the Higgs field is dominant. 

Solutions with real k represent bound states. The solution which vanishes 
at infinity is the confluent hypergeometric function: 

<t>f = x l e- kx ty(a,2 + 2l,-2kx) , a = 1 + I - - . (36) 

Equation (0) is formally the same as the Schrodinger equation for the Hy- 
drogen atom. Thus, for bound states with wavefunctions having support at 
x ^> 1, the energy eigenvalues are approximately determined by the hydrogen 
atom formula. This formula follows by demanding that the wave function 
be regular at the origin, which amounts to say that the parameter a is an 
integer < 0. Thus in this case the eigenvalues are given by 

a£ + 1 , . 

n = n + l + l = — , n = 0,1,2,... (37) 

fv 
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or 



a(n 2 - 1) ± ny/n 2 - 1\J a 2 + 1 
^ _ 2^ — 2 ' 

=* ±t/l--^ + a(l--^) , a«l . 
V n 2 n 2 

We shall see that (|38| ) is a very good approximation for all eigenstates with 
I > 1. 

Note that there are also eigenvalues of negative energies. The negative 
energy eigenstates should be interpreted in terms of the antiparticle of op- 
posite charge (thus a — > —a, see eq. (|7|) ) and positive energy. Reversing the 
sign of the energy has the same effect in the equation as reversing the sign of 
a. These states no longer correspond to the neutral system under discussion, 
but to a system of total electric charge equal to 2e (or — 2e) . They are bound 
to the dyon because the Coulomb potential is still attractive thanks to the 
Higgs contribution (i.e. a£ < but 1 + a£ > 0). Here we will discuss only 
eigenvalues of positive energies. 

Note that there are bound states even in the limit a — > 0. The origin of 
the binding force in this case is the Higgs field. 

The binding energy is represented by 

k 2 = l + 2a£-£ 2 = (£+-£)(£-£-) , 

with 



£± = a ± Vl + a 2 . 
Here represents the asymptotic mass. 
Interior solutions: For x < 1 the equations are 



[ - ^d x x 2 d x + fcll _ fcgj^ = o , (39) 

with 

k 2 n = £ 2 + — . 
The solution that is regular at r = is given by 

<l>? = c J l+h (k ia x) . (40) 
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The energy eigenvalues follow from imposing continuity of the wavefunction 
and its first derivative at x = 1. This gives the condition 



(41) 



The resulting eigenvalues (in terms of /c 2 
1. 



1 + 2a£ — £ 2 ) are given in Table 



/ 


n 


k 


k 

"'an 


^point 


1 





0.5724 


0.5953 


0.6023 


1 


1 


0.3937 


0.4048 


0.4089 


2 





0.4076 


0.4094 


0.4089 


2 


1 


0.3076 


0.3086 


0.3086 


3 





0.3086 


0.3086 


0.3086 


3 


1 


0.2476 


0.2476 


0.2476 



^num cLI"6 til 6 



Table 1: Eigenvalues for binding energy for a = 0.2. 
eigenvalues obtained by numerical integration while k an are obtained 
by the analytic method. The last column k po [ nt are the eigenvalues in 
the point-like limit, where the differential equation is formally the same 



as the hydrogen atom equation (using (38 



The numerical determination of the eigenvalues is performed by integrat- 
ing the differential equation from x = to some x = x max ^> 1, with regular 
boundary conditions at zero. The energy eigenvalues are then determined 
by the condition that the wavefunction approaches zero at x — x max (equiv- 
alently, that the wave function approaches the asymptotic form (|36"D ). The 
results are given in Table 1. Fig. 1 is a plot of the wave functions for the 
ground state and the first excited level. 

A final remark about the stability of the bound states. The theory con- 
tains also photons, i.e. the quanta of the (far from the core) unbroken U(l) 
gauge field, which are coupled to our charged scalar particles. The coupling 
is proportional to the charge operator, which has nonvanishing matrix ele- 
ments between the states a and the states S . Because S is coupled to the 
massless scalar 0Oj this implies an instability of the bound state. The domi- 
nant decay channel of the fundamental state </> a is by single photon emission, 
and subsequently conversion of S into a <po particle that can get to infinity. 
A further discussion of this is given in Section 5. 
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-0 . 2 
-0 . 4 



Figure 1: Wave functions for the ground state (solid line) and first excited state (dashed 
line) with I = 1, no = 0, 1, respectively, for a = 0.2. The dotted line is the plot of K 2 (x), 
which indicates the region of the dyon core. 



4 Scalar particles in the fundamental repre- 
sentation 

4.1 Covariant Equations 

We will work in the gauge where the charge operator is Q = ~ a 3 . Now the 
wave function ip{r) is an isospinor, that is an isospin doublet, and the charge 
operator acts on it in the standard way as a 2 x 2 matrix on an isospinor. In 
this gauge, the space components of the covariant derivatives are 

= dl + TT\ 2\[K {uiv ■ a + vin T ■ a) + vi(l - n 3 )a 3 )} . (42) 

Zr[L — n 3 ) 

Here we use the notation of Section 3.1, with U defined in eq. (|16f). 
Notice that for K = (i.e., far away from the core) one has 

A 2 = (di + ^Z ^ s) 2 
2r(l + n 3 ) 

= \d r r 2 d r + 1 (smedo(sm0d e ) + {d 9 + Ul - cos^)^) 2 ). 
r z r l sin \ 2 / 
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This is the covariant Laplace operator for a scalar particle of charge ±~e 
moving in a U(l) Dirac monopole background. 
The full equations can be written as: 



- A 2 ■*=(- ^_9y9 r + fj - ^ + ^Jl l2 )i> + j 2 M • i, , (43) 



: )■ M=±J°._ (44) 



where 

' /2= ^ o' J 2 V J ' iW= ^v-e^V 1/2 

w, = I 
Now 

Jli/2 = --^V ±1/2 V ±1/2 (46) 

with 

V±i/ 2 = sin#<% ^ id v + |(1 — cos0) , (47) 
V±i/2 =fan9d e ±id v - §(l-cos0) . (48) 



4.2 Harmonic analysis 

(± l) 

Let Z; m 2 be an eigenfunction of the operator 

Ll = J 2 ±1/2 - 1/2 + 1/4 , (49) 

In particular consider 



Explicit expressions for these eigenf unctions are given in ref. |2(J. The 
eigenvalues are now I = |, | + 1, ~ + 2, . . . Note that J. 2 .^ = + 1/2) 2 - It is 
convent to fix a phase convention by defining 



=--— gVi/aZ^, = Z + (50) 
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It is indeed easy to verify that 



Z$ = 7 , LlZ^ = Kl + l)Z$ • (51) 

\i sin U 

The angular dependence is then solved by writing 

Next, we compute • Z? = U~ l (d - |ieA )C/. Using U~ lr yU = er 3 , 
we find 

Dl • * = — (E - |egU a 3 ) 2 ^ . (52) 

Now we consider a system of vanishing total charge. To incorporate charge 
conservation, as in Section 3.2 we write the dyon charge operator as q = 
—ied x , and the wave function as 




= /+ i • (53) 



The net effect is that the cross term in the square appearing on the right 
hand side of ( |52"D (which contains the Coulomb interaction) has the same 
sign for the upper and lower components, i.e. it is an attractive potential for 
both components. 



4.3 Case of the N = 2 Hypermultiplet 

From the point of view of N = 1 supersymmetry, the N = 2 vector multiplet 
contains a iV = 1 vector multiplet and a chiral multiplet $. In this Section we 
add to the N = 2 pure Yang-Mills theory a hypermultiplet in the fundamental 
representation of SU(2). The hypermultiplet contains two chiral superfields 
Q and Q, which couple to $ by a term W = QQQ. The fundamental 
scalars will get a mass due to the coupling to the Higgs field. To simplify 
the discussion, here we will not add an independent mass term MQQ. In 
the next subsection, we will consider scalar particles with an arbitrary mass 
parameter. 

Let us consider a neutral system of dyon and scalars with charges (|, — |) 
and (— |, |). Taking into account the coupling to the Higgs field, we get the 
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following equations for 77+ (r) and r?_(r) (// being defined in (|5 
2r 2 r 2 



[--«9 r r 2 d r + — +-(^-l)] ri+ + fi T r l _ = [{E -\msG) 2 -\m 2 G 2 ]^ , (54) 



[__ 9rr ^ + _ + _( M 2_| )]?7 _ + /i_ 77+ = [(jE _l msG) 2_l m 2 G 2 ]r/ _ ^ (55) 



1 Q „ 2Q ,^! + I/..2 im_ , VK. 
2r 2 r 2 

where s = sin 6. Let us now define r] s = i] + +r/_ , r/ a = r/ + — 77_ . They satisfy 
the decoupled equations 

'2 



1 K 2 1 (//<■" 

[-^ 2 ^+^ + ^(^ 2 -|)]^+72-^ = [(^->^) 2 -> 2 G> s , (56) 



[-^ 2 ^ + ^ + 3(^-DK-^ a = [(E-|m S G) 2 -> 2 G> a . (57) 



-5 r r 2 «9 r + — + -(//-§)K-^ ■ 

Note that the equations for r? s and 77 a are formally the same under the ex- 
change /i — > — /i, so they can be investigated on the same footing. In what 
follows we consider the 7] s equation (|5"6), wherefrom we obtain the solutions 
for r7 a by flipping the sign of /i. 

In terms of x = (m cos 6) r, the equation fl56|) for 77 s reads 

r 1 o 20 -^ 2 (a;) 1 , 2 1N , /iif(x) 

[ - - 2 d x x 2 d x + + - 2 ^ 2 - \)K + ^^vs 

= [£ 2 - a£G(x) - \G\x)} Vs , (58) 

£ = — , a = tan 6 . 

mcos O 



At x ^> 1, the differential equation (|58"D takes the form 

-^d x x 2 d x r] s + -^-^ + k 2 r/ s = , (59) 



with 



± + a£-£ 2 . 



4 

The solution which vanishes at infinity is the confluent hypergeometric func- 
tion: 

Vs = x l e- kx ^(a,2 + 2l,-2kx) , (60) 
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1 + / 



2a£ + l 
4k 



The approximate energy eigenvalues are then determined by the condition 

2a£ + 1 



n + I + 1 



4j\ + a£-£ 2 



n 



0,1,2,... 



(61) 



Some values are given in Table 2. Note that energy eigenstates obtained in 
this approximation are the same for rj s and r) a , since at x 3> 1 r/ a satisfies the 
same eq. (0). This point-like approximation is already good for states with 
1 = 3/2 and becomes better for states with high angular momentum, which 
are farther from the core. 



1 


n 


^imm (s) 


^num ('V 


h ■ 

""point 


1/2 





0.1921 


0.2245 


0.2045 


1/2 


1 


0.1188 


0.1310 


0.1238 


3/2 





0.1236 


0.1241 


0.1238 


3/2 


1 


0.08848 


0.08881 


0.08863 


5/2 





0.08863 


0.08863 


0.08863 


5/2 


1 


0.069005 


0.069005 


0.069005 



Table 2: Eigenvalues for binding energies /c nU m(s,a) for r] s and r) & using 
a = 0.2. fcpoint are the approximate eigenvalues given by the analytic 
expression (|61|). 



At x -C 1, the equation is 
1 



.i- 



d x x d x r] s + 



with 



£ 2 + 



6 



x- 



(62) 



l + h 



The solutions are Bessel functions. 

In the case of ^ a , the "centrifugal barrier" for x <^ 1 is . As a result, 
there is no barrier for the ground state with fi = 1, and the corresponding 
wave function does not vanish at the origin (see fig. 2). 

The numerical determination of the eigenvalues is performed as in Section 
3 by integrating the differential equation from x = to some x = x max ^> 



18 



0.8 



0.4 



0.2 



0.6 



1 




5 



10 



15 



20 



25 



Figure 2: Wave function for the ground state r] a (solid line) and for % (dashed line) with 
I = 1/2, n = (we use a = 0.2). The dotted line is a plot of K 2 , indicating the region of 
the dyon core. 

1, with regular boundary conditions at zero. The energy eigenvalues are 
then determined by the normalizability requirement that the wavefunction 
approaches zero at x — x max . The results are given in Table 2. 

4.4 Massive scalars in N < 1 without coupling to Higgs 

In the systems of Section 3.4 and 4.3, for small a the Coulomb attraction is 
dominated by the Higgs field. This interaction is dictated by iV = 2 super- 
symmetry. As a result, the bound states are not very sensitive to the value 
of an a < 1. It is of interest to investigate systems with less supersymmetry, 
to see the effect of changing the values of the coupling. Here we will consider 
a scalar field in the fundamental representation with an arbitrary mass term 
and no coupling to the Higgs field. 

Let us consider a neutral system of dyon and these scalars with charges 
(|, — |) and (— |, \)- In a similar way as in Section 4.3, we find the following 
equations for 77+ (r) and ?7_(r) : 



i-^dr + ^ + ^-^ + ^V- = [{E-\msGf-M^ + , (63) 
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1 K 2 1 11 K 

l-^rr% + — 2 +^-^)v- + ^V + = [(E — \msG) 2 — M 2 ]rj- , (64) 

where we have introduced a mass M. Thus, again the presence of r/ + turns on 
the other component (and viceversa), so the bound state is an entangled 
quantum state. The equations are decoupled in terms of i] s = i] + + r\- , 
r) a — r) + — ?7_, which satisfy 

1 K 2 1 a K 

[ " ~2 d rr 2 dr + ^ + ^ - + ^Vs = l(E - \msGf - M 2 ]rj s , (65) 



1 K 2 1 a K 

[ " ~M d r + ^2+-2^- ~ = K E ~ I™ )" ~ • (66) 

Introducing the rescaled radial coordinate x — (m cos 6) r, the equation for 
rj s becomes 



r 1 _ 2 „ K 2 (x) 1 2 , uK(x) 

[ - -d x X 2 d x + -±± + - 2 (^- \)\Tk + 

= [(£- §aG(s)) 2 - M 2 ^ , (67) 



£ = — , M 



m cos 6 m cos 

In the region x >> 1, this differential equation reduces to 



- ^d x ^d xVB + (1(^ 2 - I - I a 2 ) - ^ + k 2 ) Vs = , (68) 

OC \ X X / 



with 



k = ^M 2 - {£ 



-la) 2 



The solution which is normalizable at infinity is the following confluent hy- 
pergeometric function: 

Vs x b e- kx ty(a, 1 + c, -2kx) , (69) 
6=i(_l + c ) > a= i(i + c )_g, C=v / V _!_ a 2. 

The approximate energy eigenvalues in this point-like limit are then deter- 
mined by the formula 



= 1(1 + yV -I -a 2 ) aS , n = 0, 1, 2, ... (70) 

z v v o \A2 C 1^A2 



n 



2 _ (5 _ i a)2 
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Some eigenvalues are given in Table 3. Note that the binding energies are 
very small compared to the asymptotic mass e x = M + fa. This can be 
understood as follows. From eq. ( [7C| ) we see that 



AS = 8^-8 

This is of 0(a 2 ) if M 
li = 1, no = and M. - 



a 2 {2M + af 
8M(2n + l + c) 2 

-- 0(1), and of 0{a 3 ) if M 
a one has A8 = 0.15 a 3 . 



8 n 



M + \a 



0(a) 



(71) 

For example, for 



1 


n 


^num(s) 


^mim ( 9>) 


h ■ 

"'point 


1/2 





0.03626 


0.03728 


0.03666 


1/2 


1 


0.02104 


0.02140 


0.02117 


3/2 





0.02053 


0.02053 


0.02053 


3/2 


1 


0.014553 


0.014555 


0.014545 


5/2 





0.0144460 


0.01460 


0.014460 


5/2 


1 


0.011217 


0.011217 


0.011217 



Table 3: Eigenvalues for binding energies /c num (s,a) for rj s and r\. A using 
0.2 and M. = 0.5. k point are the approximate eigenvalues given by 



a 



the analytic expression ([70|). 



In the region x C 1, the equation becomes 

-^^+(^±^-^ = 0, (72) 

with 

6 

The solutions are Bessel functions. 

The numerical calculation of the eigenvalues is performed as in the pre- 
vious cases by integrating the differential equation from x = up to some 
x = x max 3> 1, with regular boundary conditions at zero. The energy eigen- 
values then follow from the requirement that the wavefunction approaches 
zero at x = x max . Some results are given in Table 3. 

As in the case of the previous subsection, for the states of 7/ a with \i = 1 
(i.e. I = 5)) the potential has no centrifugal barrier at x = 0, and the 
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corresponding wave functions do not vanish at the origin. The equation in 
the interior 2; < 1 becomes 

- \d x x 2 d xVs - kl Vs = , (73) 
ar 

with 

i& = -M 2 + e 2 . 

Although there is no centrifugal barrier near x — 0, it is worth noting that 
the particle is not concentrated at the core. Its extension grows like 

Ax 



AS Ma 2 ' 

which is much greater than one for small coupling a. 

The wave functions for the lowest energy states of ?7 S and 7] & are shown in 
fig. 3, for the case a = 0.2, M. = 0.5. Comparing to the case of fig. 2, where 
the Coulomb interaction is dominated by the coupling to the Higgs field, we 
see that the wave function is more extended. This is expected, in view of the 
above observation that the extension is greater for smaller couplings. 

The lowest energy state for i] s should be unstable because the 7/ s entangled 
state can make a transition to a r? a entangled state by emission of a photon 
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(see also the discussion in Section 3). However, the lowest energy bound 
state for i] a (with I — |, no — 0) should be stable, since there is no possibility 
of decay into a state of lower energy. The same conclusion applies for the 
bound states of Section 4.3 . 

5 Summary and Discussion 

The main problem to which we addressed our study was the nature of the 
possible bound states of a scalar particle around a dyon, in systems which 
are globally electrically neutral. Formally, this amounts to considering the 
scalar field (different from the Higgs) as a fluctuation in the background of 
the dyon solution, and treating it in the second quantization formalism. Since 
we have taken as a paradigm of the dyon the Prasad- Sommerfeld solution of 
N = 2 supersymmetric SU (2) YM theory, we first considered scalar fields in 
N = 2 supermultiplets. In these cases, the scalar particle has an attractive 
interaction with the long range Higgs field, whose strength is fixed by the 
dyon solution, besides the usual Coulomb force (which is proportional to the 
square of the electric charge, and we naturally take it to be quite weaker 
than the other). In Section 4.4, we considered a case where the scalar field 
does not interact with the Higgs, which can be the case in systems with 
less supersymmetry. Here the bound state is solely due to the Coulomb 
attraction. The main difference with respect to the previous cases is that 
bound states are, in this case, much more extended far from the monopole 
core. This is expected, since the Coulomb interaction is weaker. 

It should be noted that the bound states are rather larger than the 
monopole core even for the cases where the Higgs attraction is dominant. 
Therefore they are expected to be robust against back-reaction effects of the 
classical dyon solution for sufficiently small coupling constant (apart from 
the effect of charge conservation which we have already included). 

In the case of the scalar field in the adjoint representation, we have con- 
sidered it to be a member of the N = 2 vector multiplet. Therefore it 
is coupled to the Higgs and it takes mass through the Higgs expectation 
value. The attractive force is dominated by the Higgs field and there are 
bound states also in the limit of vanishing charge. The bound states corre- 
spond to a certain linear combination |dyon + ) (scalar _) — |dyon_) |scalar + ), 
whereas the other linear combination |dyon + ) (scalar _) + |dyon_) (scalar + ) 
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mixes with the state |monopole ) |scalaro), in which both the particle and 
the monopole are uncharged. This mixing cannot form bound states because 
the neutral particle is massless and therefore it can have an arbitrarily low 
energy and still escape to infinity. Since some perturbation not included in 
our analysis, say the e.m. quanta radiation, could cause a transition from 
|dyon + ) |scalar_) — |dyon_) |scalar + ) to |dyon + ) | scalar.) + |dyon_) |scalar + ), 
we conclude that these bound states are unstable. 

In the case of the scalar field in the fundamental representation, we have 
considered two cases. In the first case, the scalar is in a iV = 2 hypermultiplet 
and is coupled to the Higgs field by maintaining the N = 2 supersymmetry 
of the Lagrangian. We have assumed that its mass is completely due to the 
Higgs expectation value, in order to compare with the previous case. The 
resulting bound states are again essentially due to the Higgs attraction. In the 
second case, we have explicitly studied a SUSY breaking scenario, in which 
the coupling to Higgs is absent and there is an arbitrary mass parameter. 
The qualitative features of the two cases are rather similar. 

Since in the fundamental representation there are no particles of zero 
charge, both combinations |dyon + ) (scalar _) ± |dyon_) |scalar + ) admit bound 
states (here we consider dyons with the same ±e/2 charges as the particles in 
the fundamental representation, in order to have a globally neutral system). 
The lowest energy level is of the form |dyon + ) | scalar _) — |dyon_) |scalar + ), 
with the angular momentum taking the minimum value I — \. In this case 
the "interior centrifugal barrier" is weak and the wave function is nonzero at 
the origin, where typically it takes its maximum value (an exception occurs 
in the second case with no coupling to Higgs, when the charge is very small: 
in such case, we found that the maximum of the wave function appears at 
some finite r). The extension of the wave function is in both cases larger than 
the monopole core size. The orbiting particle can thus be always considered 
to be mostly outside the core, with a wave function similar to that of the 
Hydrogen atom. The main role of the monopole core is to produce a full 
quantum entanglement, despite its small effect on the r dependence of the 
wave function. 

Thus a bona-fide bound state picture emerges in both cases. We con- 
clude that we get stable, quantum entangled, bound states in the case of the 
fundamental representation. 

While generic quantum entanglement of charges may occur in various 
physical situations (for instance in the case of the decay of a neutral particle 



24 



in two charged components), the remarkable phenomenon found here is that 
stable bound states of a light particle and a heavy dyon can only exist in 
an entangled state involving the system with opposite charge. The bound 
particle and the dyon can be very far from each other, and they are not charge 
eigenstates (although the total charge is zero). Although this setting does 
not seem to be in conflict with any physical law, it is nevertheless a curious 
effect and one could amuse her- (him-) self by imagining would be paradoxes 
in the chemistry of such an "atom" . 

Note that a measurement of the charge of the orbiting particle should 
make the wavefunction precipitate into a defined charge eigenstate (say, from 
|dyon + ) |scalar_) — |dyon„) |scalar + ) to |dyon + ) |scalar_)). Since such bound 
state is not possible, one concludes that any measurement of the charge 
should require an energy above the binding energy, so that the final scalar 
particle can escape to infinity. 

Similar bound states should exist for fermions and vector particles i.e. 
they should give rise to a state consisting of a mixture of fermions or W + 
and W~ gauge bosons entangled with dyons. This may be more easily derived 
by embedding the N = 2 supersymmetric model in iV = 4 super Yang-Mills 
theory, so that the solutions of the wave equation for fermion and gauge boson 
fluctuations are connected by an unbroken supersymmetry to the solutions 
for scalar fields computed here. 
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7 Appendix: Harmonic Analysis for adjoint 
scalars 

In this Appendix we perform the harmonic analysis of the coupled equations 
for the case of scalar particles in the adjoint representation (the harmonic 
analysis in the case of a pointlike monopole has been done in ref. ||19|). 
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Using eqs. fllTD - (0), we obtain the following system of equations in 
spherical coordinates: 



K , . _ _ „ . K 



OC- 



sinOde - id^Fo - a + . 2 (sinflgg + id v )F , (74) 



r z sin r z sin # 



p ±iip 1 1 1 - 

^•[d T F ±] = ^—[-dyd r --(-i + K*)--ji) h± 

suit' 

where ^ 

L 2 = — r~T7:(sm9de + idn,)(sm6d e - mL) , (76) 
sin 

is the standard angular momentum and 

Jl = — ^V + V + , J 2 = — 1-V_V_ , (77) 
sin sin 

is a covariant angular momentum. Here 

V± = sin Ode =F id v + (1 — cos 0) , 

V± = sin Ode ± id v — (1 — cos#) , 
[V±,V±] = -2sin 2 # . 
The full equations of motion can now be rewritten as 

( - - 2 d r r%. + % + - 2 LAf - -4^i^+ h + + e-*^-h-) = E 2 F , 
V r z r z r z J r A sin 

(78) 

[ - ^d r r 2 d r + £(-1 + K 2 ) + ±j 2 }h + + ^(sm9d e - id v )F 

= [(E - eVq) 2 - e 2 a 2 G 2 }h + , (79) 



[ - ^d r r 2 d r + £(-1 + K 2 ) + ±J 2 }h_ + "^(sinedo + id 9 )F, 

= [{E + eVq) 2 - e 2 a 2 G 2 ]h_ . (80) 
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The angular dependence is solved by setting 

F o (f) = <f )o (r)Y lm (0,v) , (81) 

h+(f) = y Mr) z+je, <p) , M*0 = f <^ r ) <p) > ( 82 ) 

'0 '0 



Z = V«0 + 1) , (83) 
where (p) are the spherical harmonics, i.e. 

L 2 Y lm = 1(1 + l)Y lm , (84) 
with / = 0, 1, 2 . . ., and Z lrn are eigenfunctions of the operator J 2 , 

J 2 Z, m = 1(1 + l)Z lm . (85) 
Explicitly, they are given by 

ZL = -^(sm ^ e =F id v )Y lm . (86) 
sm (7 

Using ([86"D , one can check that ( jS5|) is indeed satisfied. 

Here / = 1,2,.... Note that there is no eigenfunction of J 2 for I = 0, 
since there are no normalizable solutions of V±Zq = 0. 



Using eqs. (|7q)-(|82[) and the properties: 

itp _ -iip 

-V+Z+ = -1(1 + l)Y lm , -t-^V-Z,- = -Z(Z + l)^ m 



sin# T lm v ' ' sinfl 
we find eqs. ©-(H). 

The functions Z^ (|86| ) have a simple form in terms of 9, For example, 
for 1 = 1 we have the following solutions 

to = : Fo(f) = <p (r)cos9 , h±(f) = —= <fi±(r)e Tlip sin 9 

v2 

m = 1 : Fo(r) = 4>o(r)e tip sin9 , h + (r) = —=<j) + (f)(\ + cos^) , 

V2 

/i_(r) = -^</>_( r )e 2 ^(-l +cosfl) , 
v2 

to = -1 : F (f) = o (r)e-^sin^ , /i+(f) = -^</>+(r) e - 2 ^(-l + cosfl) , 

v2 

/i_(f) = — 0_( r )(l + cos6>) . 
v2 

Similarly, one can write down expressions for higher I. 
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